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We review recent work on the Bose-Einstein condensation of photons in a dye mi- 
crocavity environment. Other than for material particles, as e.g. cold atomic Bose 
gases, photons usually do not condense at low temperatures. For Planck's blackbody 
radiation, the most ubiquitous Bose gas, photon number and temperature are not 
independently tunable and at low temperatures the photons simply disappear in the 
system's walls, instead of massively occupying the cavity ground mode. In the here 
described approach, this obstacle is overcome by a fluorescence-induced thermaliza- 
tion mechanism in a dye-filled microcavity. Experimentally, both the thermalization 
of the photon gas and, at high photon densities, Bose-Einstein condensation has been 
observed. This article describes the thermalization mechanism of the photon gas in 
detail and summarizes so far performed experimental work. 
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I. INTRODUCTION 

When a gas of particles with given density is cooled to such low temperatures that the 
associated de Broglie wavepackets spatially overlap, quantum statistical effects come into play. 
Specifically, for a gas of particles with integer spin (bosons) , Bose-Einstein condensation into 
the ground state sets in above a critical phase space density. For dilute atomic gases, this 
effect has first been observed in 1995 by means of laser and subsequent evaporative cooling 
of alkali atoms il-3|, see also the article of F. Chevy and J. Dalibard in this volume. The 
signature of Bose-Einstein condensation has also been observed for several solid-state quasi- 
particles, as exciton-polaritons and magnons, see e.g. the contributions by Y. Yamamoto and 
by V. Demokritov and A. Slavin, respectively, in this volume. 

Photons, the quantized particles of light also are bosons, but usually show no Bose-Einstein 
condensation. In a blackbody radiator, the chemical potential of photons vanishes, i.e. the 
average particle number does not follow a given conservation law, but adjusts itself to the 
available thermal energy ^. This is the essence of the Stefan-Boltzmann law, linking the total 
radiation energy U to the fourth power of the temperature, U oc T^. At low temperatures, 
the photon number simply decreases and no macroscopic occupation of the cavity ground 
state occurs. Thus, a necessary precondition for a Bose-Einstein condensation of photons is 
to find a thermalization process that allows for an independent adjustment of both photon 
number and temperature. The lack of such a mechanism has long prevented the realization 
of light sources that are capable of generating single mode light, without the necessity to be 
driven out of thermal equilibrium, as e.g. for a laser. Note that in a laser both the state of 
the light field and that of the active medium are far removed from thermal equilibrium j^. 
To some extent, lasing is even a prime example for a non-equilibrium process, as only the 
absence of thermal equilibrium allows for inversion and optical gain. Early theoretical work 
has proposed to reach Bose-Einstein condensation of photons by Compton scattering off a 



thermal electron gas Later, Chiao et al. proposed a two-dimensional photon quantum fluid 
in a nonlinear Fabry-Perot resonator, where thermalization was sought from photon-photon 

g processes in atomic physics 

BEC experiments, though the limited non-linearity has so far prevented a thermalization of 
the photon gas |10| . In other work, the demonstration of (quasi-)equilibrium Bose-Einstein 
condensation of exciton-polaritons, mixed states of matter and light, has been reported 



13| . Here interparticle collisions of the excitons, i.e. the material parts of the polaritons, act 
as a therma 
observed |l4. Il5|. 



ization mechanism. In other experiments, superfluidity of polaritons has been 



In recent experiments of our group, photon Bose-Einstein condensation is achieved in a 



dye-solution fllled optical microresonator 



Thermalization of the photon gas with the 



dye is achieved by repeated absorption emission cycles. For such systems it is known that fre- 
quent collisions (~ 10 fs timescale) between solvent and dye molecules causes rapid transverse 
decoherence at room temperature, so that the condition of strong light-matter coupling is 
not met H Q. The distance between the two spherically curved resonator mirrors is in the 
micrometer regime, which causes a large frequency spacing between the longitudinal resonator 
modes. The latter is of order of the emission width of the dye molecules. In combination 
with an intracavity modiflcation of the spontaneous emission, preferring the emission to small 
volume modes (low transversal excitation), a regime is reached, where to good approximation 
the resonator is populated only by photons of a single longitudinal mode number, see Fig. 
[1^. The longitudinal modal quantum number is frozen out and the photon gas effectively 
becomes two-dimensional. As is indicated in Fig. [TJ), the photon dispersion relation acquires 
quadratic, i.e. particle-like character, with the frequency of the transverse TEMqo mode act- 
ing as a low-frequency cutoff frequency. Furthermore, a harmonic trapping potential for the 
photon gas is induced by the mirrors curvature. Thermal equilibrium of the photon gas with 
its environment (at room temperature) is achieved as the photons are absorbed and emitted 
by the dye molecules repeatedly. The photon frequencies will accumulate within a spectral 
range of order ~ ksT/h above the low frequency cutoff. Other than in a blackbody radiator, 
the thermalization process allows for an independent adjustment of temperature and photon 
number. This becomes clear by noting that the energy of fluorescence photons, which takes 
values above the cutoff frequency of the resonator in energy units (/itUcutofi — 2.1 eV), is far 
above thermal energy, k-^T ^ 1/40 eV. Thus, purely thermal excitation of (optical) photons is 
negligible. Instead, the number of photons will be determined by the strength of the optical 
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Figure 1: (a) Scheme of optical resonator (top) and the cavity modes (bottom) for the case of a 
spacing between resonator mirrors of half an optical wavelength of the lowest optical mode. The 
resonator is filled with dye solution. The photon gas thermalizes to the temperature of the dye 
solution, (b) Photon dispersion in resonator (solid line) and the dispersion of a free photon (dashed 
line). 



pumping. One can show that the photon gas confined in the resonator is formally equiv- 
alent to a harmonically trapped two-dimensional gas of massive bosons with effective mass 
meff = hWcntoRnQ/c^ , whcrc no denotes the refractive index of the medium and c the vacuum 
speed of light. For such a system it is well known that a Bose- Einstein condensate exists at 



a finite temperature 2l|, |22|. In recent experiments, we have both observed thermalization of 
;he photon gas in the dye-filled microcavity system 17| as well as Bose-Einstein condensation 



16|. The so far observed properties of the photon Bose-Einstein condensates in many re- 



spects resemble that of atomic gases, while the approximately 10 orders of magnitude smaller 
effective photon mass allows for transition temperatures in the room temperature regime. 

Bose-Einstein condensation and superfiuidity are two closely related phenomena 2^. 
While the former is connected to equilibrium properties, the latter deals with transport prop- 
erties. Bose-Einstein condensation is in principle possible with an ideal gas, while the presence 
of superfiuidity requires interparticle interactions. In future, it remains to be experimentally 
verified whether the photon Bose-Einstein condensate also exhibits superfiuidity. 

In the following, section II describes the fiuorescence induced thermalization mechanism 
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of the photon gas and section III the statistical theory of the trapped photon gas. Further, 
section IV reviews experiments on the thermalization process and section V corresponding 
results on Bose-Einstein condensation. Finally, section VI gives conclusions. 



II. FLUORESCENCE INDUCED THERMALIZATION 
A. Kennard-Stepanov theory of dye spectra 

Early experimental work has shown that spectra of dye molecules in liquid solution show 



several universal properties j24j-|27j. These are for example the mirror rule, which states that 
the spectral absorption profile is a mirror image of the fluorescence profile; the Stokes rule, 
stating that the spectral centroid of fluorescence occurs at a higher wavelength than that of 
absorption; and Kasha's rule [2b|, which expresses that the fluorescence does not depend on 
the wavelength of the exciting light. These common properties, which to good accuracy are 
fulfilled in many dye species, suggest that absorption and fluorescence in such systems follow 
a general mechanism. Most of these properties can be understood as being the consequence 
of a collisionally induced thermalization mechanism 2J, |25|. For a corresponding model. 



consider an idealized dye molecule with an electronic ground state 5*0 and an electronically 
excited state Si, each of which are subject to additional rovibronic level splitting, as shown 
in Fig. [2^. The vibrational and rotational state of a dye molecule in liquid solution is 
permanently altered by collisions with solvent molecules, e.g. on a femtosecond timescale at 
room temperature. The frequent collisions lead to a thermalization of the rovibrational state 
within a sub-picosecond timescale, which is much faster than the electronic lifetime of the 
excited state, being typically of the order of nanoseconds. When a photon is absorbed, the 
dye molecule is likely to be transfered to a highly rovibronically excited substate of the Si 
manifold, but the excessive rovibronic energy will be quickly dissipated into the solvent bath. 
The fluorescence photon will be emitted from a dye state that is in thermal equilibrium with 
the solvent bath, with typically lower rovibrational quantum number. This thermalization 
process, which occurs both in the electronically excited and in the ground state level, explains 
why fluorescence can be dissipative (Stokes shift) and why there typically is no correlation 
between the wavelength of the absorbed and the emitted photon (Kasha's rule). 

The thermalization of the rovibronic degrees of freedom has another, though closely related, 
consequence: the Einstein coefficients of absorption and emission at a certain photon energy 
hw are connected by the Boltzmann factor of that energy. This relation is known as the 
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Kennard-Stepanov law, and can be written in the form 



ft("-"o) 



(1) 



where Bi2{uj) and i?2i(w) are the Einstein coefficients of absorption and stimulated emission 
respectively, Uq is the frequency of the zero-phonon line of the dye, and Wi^-^ are statistical 
weights related to the rovibronic density of states, which will be defined subsequently in this 
text. We note that the Kennard-Stepanov relation can also be stated in terms of A{uj) / Bui^u)), 
where A{u) is the Einstein coefficient for spontaneous emission - with the difference to the 
above definition essentially being the density of states. This relation has been discovered in 



the beginning of the last century 
historical outline can be found in 
be found in the literature 31 



29 



30|, and has been 'rediscovered' several times. A short 



31] . Both theoretical and experimental investigations can 



In the following, we give a short derivation of the Kennard-Stepanov relation. As described 
above, we model the dye molecule by an electronic two-level system with levels 5*0 and Si, 
also denoted by J, and f, each of which is subject to additional rovibronic level splitting. It is 
important to note that the Einstein coefficients of such a medium at a given frequency u are 
an average over all pairs of individual rovibronic substates with a G 5*0, P E Si that 

match the transition frequency: 



ea + fiw = tvjjQ + ep . (2) 

The latter equation expresses energy conservation, see also the Jablonski diagram of Fig. 
We assume that the population of rovibronic states in both lower and upper electronic states 
is fully thermalized from frequent collisions with solvent molecules. This assumption will be 
valid as long as the radiative lifetime of the electronically excited state 5*1 remains clearly 
longer than the thermalization time. We thus expect that the corresponding substates within 
the lower and upper electronic manifolds will be occupied with a propability given by the 
Boltzmann factors 

Pa = G ''b'^ /wi and Pi3 = ^ '-B^/w^ , (3) 
with the normalization factors 

W]^ = ^ and = ^ . (4) 

aeSo i3eSi 

From eqn ([2]) and ([3]) one immediately obtains 



7 



(a) 



Si 







s 


— 1 






fico 






-J 


' 


a 



lea 



(b) 




ground state optical 
molecule j photon 

(So) 




excited 
molecule \ 

(Si) 



Figure 2: (a) Jablonski diagram of a dye molecule with two electrouic levels. Both lower and upper 
electronic state, Sq and Si, are split into a rovibronic substructure, (b) The process of absorption 
and fluorescence of a photon can be seen as a photochemical reaction. 



P/3 = — e '=bt p„ . (5) 



We can now write 



where we have introduced the Einstein coefficients B{a — )■ /3) for transitions between the 
corresponding individual rovibronic states. Applying B{a — (3) = B{(i — )■ a) and using in 
addition ([5]), we finally obtain 

= ^e-"^ , ' (7) 

which is the Kennard-Stepanov law of eqn ([1]). As described above, the Kennard-Stepanov 
relation is experimentally well established. However, many dye species show smaller or even 
larger deviations. For a discussion of these issues, see section IIIEI 



B. Chemical equilibrium between photons and molecules 

Molecules in the electronic ground state can be transferred to the upper electronic level 
by the absorption of a photon, while photon emission leads to de-excitation. These processes 
can be seen as a photochemical reaction of the type 



(8) 
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see also the illustration of Fig. I^b- Here 7 stands for a photon, J, for a molecule in the electronic 
ground state, and for a molecule in the electronically excited state. In chemical equilibrium 
the corresponding chemical potentials, with fi^ (fi^) denoting the chemical potential of a 
ground (excited) state molecule and fij indicating the photon chemical potential, satisfy the 
relation 

fi^ + fi^ = fi^. (9) 
If we express this relation in terms of the photon fugacity z, the equation reads 

Next, let the partition function of a dye molecule be 

where w^^-t- are the statistical weights of eqn Further, ujq is the frequency of the zero- 
phonon line of the dye. We can then identify 

wte /jt = ^ ^ (12) 

w^e^-^-r /-^ = p ' ^^^^ 

as the probability of finding a molecule in the excited (ground) state, where (p|) is the 
spatial density of the excited (ground) state dye molecules respectively, and p is the total dye 
molecular density. Using eqn (fTOj) . we find that the photon chemical potential is determined 
by the relation 

z = e^'^T = -t e^'BT . (14) 
^^t Pi 

In chemical equilibrium the photon chemical potential is thus determined by the excitation 
ratio p^/pii which corresponds to the relative number of dye molecules in the electronically 
excited state llSl. 13611 . 



C. Thermal equilibrium and Markov-processes 



To begin the discussion of the thermalization process, let us characterize thermal equi- 
librium in terms of a Markov process. In general, the time evolution of a physical system 



coupled to a heat bath can be regarded as a random walk in configuration space 



37 



38|. Such 



a random walk is fully characterized by the present state of the system and fixed transition 
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rates between configurations (process 'without memory'). Not all sets of transition rates are 
physically meaningful. A necessary requirement is that for a sufficiently long time evolution 
of the system each configuration has to occur with the statistical weight of its Boltzmann 
factor. In the following, we derive a condition from which one can easily decide if a given set 
of transition rates leads to thermal equilibrium. In the next section this will then be applied 
to the fiuorescence induced thermalization process of a photon gas. 

Let K denote a certain state of a physical system (for a photon gas this will be given by 
a set of mode occupation numbers). Furthermore, pxit) denotes the probability to find state 
K at time t. The temporal evolution of pxit) is given by the master equation 

PK{t + 1) - PK{t) = J]px'(^) RiK' -^K)~ J]pi^(t) R{K ^ K') , (15) 

K' K' 

where the coefficients R{K — )■ K') denote the transition rates between configurations. We are 
interested in transition rates that asymptotically bring the system into thermal equilibrium, 
i.e. which yield pxit — oo) = exp(— i^/^/^B^)/^) with Z as the partition function. The 
asymptotic master equation for this case is given by 

= ^ expi-EK'/k^T) R{K' ^ K) -J^ exp{-EK/kBT) R{K K') . (16) 

K' K' 

Equation (fT6|) has many solutions. One solution stands out by showing no net probability 
fiow between two given states (detailed balance): 

= exp(-^K'/fcBT) R{K' ^K)- exp^-EK/k^T) R{K K') yK,K' , (17) 



or 



R{K ^ K') _ g^p(_^^/^^^) y ^/ 



R{K' -> K) 

with an energy difference AE = Ek' — Ek- If the transition rates R{K — )■ K') fulfill eqn 
(|T8|) . this gives a sufficient condition that the Markov process drives the system into thermal 
equilibrium. 



D. Light-matter thermalization process 

To begin with, consider an experimental scenario as depicted in Fig. [3l consisting of a 
macroscopic box with refiecting walls that is filled with a dye solution. Light, emitted e.g. 
by a spectrally narrow optical source, is irradiated at point A. What is the spectrum of the 
radiation leaving the dye container at point B? If the optical density of the dye is sufficiently 
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Figure 3: Multiple absorption-fluorescence cycles in a macroscopic dye-filled photon box. Asymp- 
totically, the photons entering the box at point A would be transferred to blackbody radiation at 
the temperature of the environment. However, without a low-frequency cutoff the thermalization 
process remains incomplete, as photons leave the absorption bandwidth of the dye. 



large, multiple absorption-emission fluorescence cycles will occur. One then experimentally 
finds that the spectrum at point B is red shifted with respect to that obtained for a single 



fluorescence event, 
inner filter effect 



24 



n the context of fluorescence spectroscopy, this effect is known as the 
2^. Interestingly, this red shift of the fluorescence can be understood 
in terms of a partial thermalization process. By fluorescence and reabsorption, the photons 
exchange energy with the dye solution. This energy exchange establishes a thermal con- 
tact between the light and the dye medium. In principle, one would expect that the state 
of the light field asymptotically relaxes towards a thermal state at the temperature of the 
dye solution, i.e. after many complete cycles the spectrum of blackbody radiation at room 
temperature would be obtained. However, owing to the succesive red shift, along with the 
limited spectral bandwith of the dye, the reabsorption probability rapidly decreases, and the 
dye solution gradually becomes transparent. The thermal contact soon breaks down and the 
process remains incomplete. Such an incomplete thermalization occurs e.g. in some solar 
light concentrators jsQ]. 

Regarding photon number conservation, it on the other hand is not desirable to transfer the 
irradiated light into blackbody radiation. As soon as the photon energy becomes comparable 
to thermal energy, one would lose the ability to experimentally tune the photon number. In 
equilibrium, the photon number would then readjust following the available thermal energy, 
independently of the optical pumping power irradiated at point A (provided that the pumping 
does not significantly heat up the dye solution). However, in a microcavity environment, as it 
is used in our experiment, both a full thermalization and an experimentally tunable photon 
number can be achieved. This becomes clear by noting that the low-frequency cutoff, imposed 
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by the resonator, prevents a succesive red-shift of the photon gas. The thermahzation by 
repeated absorption-emission cycles can now proceed until the photon gas is fully thermalized, 
i.e. thermal contact to the dye solution is not lost. Owing to the low-frequency cutoff, the 
photon energy can remain far above the thermal energy. With ^cutoff — 2.1eV ^ 1/40 eV ~ 
this limit is fulfilled in our experiment. In this situation, no photons are created or 
destroyed on average. The thermahzation process here conserves the photon number, with 
the latter being determined by the strength of the optical pumping. 

In the remainder of this section, we present a more rigorous treatment of the light-matter 
thermahzation process is considered as a random walk in the 

configuration space of all allowed light field states. Here a state K is given by the cavity mode 
occupation numbers K = (n^,nf^,n^, . . .). The mode occupation numbers are permanently 
altered by photon absorption and emission processes. In first order perturbation theory, the 
rates (per volume) for absorption and emission of one photon in mode i at cavity position r, 
denoted by /2^'*(r) and -R^'*(r), have the form 

i?f/(r) = B,,{co,)u\r)p^nf (19) 
<'Xr) = B,,{uj,)u\r)p^{nf + 1) (20) 

where u*(r) is the spectral energy density of one photon in mode i. We assume that the 
number of dye molecules is sufficiently large, such that photon absorption and emission leave 
the densities of ground state (p^) and excited state (p^) molecules unchanged, i.e. p^^^ can 
be treated as fixed parameters. This assumption corresponds a grandcanonical limit, see also 
section IIII CI 

Suppose that a state K' emerges from state K by the absorption of a photon in mode 
i, with nf = nf — 1, and accordingly, that K emerges from K' by an emission process 
into this mode. The corresponding rates are i?^'*(r) and i?^'*(r), and with eqn 0191) and 
(120]) . their ratio is given by /2^'*(r)//2^ '*(r) = Bi2{uJi) Pi/ B2i{uJi) p^. In section III CI it was 
shown that the random walk given by the transition rates of eqn (fT9|) and ( 120|) will lead to 
thermal equilibrium, if the ratio f2^'*(r)/i?^ '*(r) is given by the Boltzmann factor of the 



energy difference between K and K' 

Bl2{uJi) pi 

B2i{uJi) Pt 



(21) 



see also 
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40 



4l|. If we now apply the Kennard-Stepanov relation, eqn ([T]), and assume 



chemical equilibrium, eqn (1141) . the detailed balance condition of eqn ( 12T]) is indeed verified. 
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This means that multiple absorption-emission cycles drive the photon gas into thermal equi- 
librium with the dye solution at temperature T, and with a photon chemical potential yU^ 
determined by the molecular excitation ratio. The average occupation number Ui of mode i 
can be determined by balancing the average absorption and emission rates at a given cav- 
ity position. As expected, this gives a Bose-Einstein distribution for the average occupation 
number, with nj = (exp [{fiWi — / k-QT] — l) ^ . 



E. Spectral temperature 

From the Kennard-Stepanov theory of dye spectra, as presented in section [Tl A I one expects 
that the spectra of absorption and emission are interlinked by a Boltzmann factor. In general, 
this law is empirically well confirmed. However, it is known that many real dyes show more or 
less pronounced deviations. Before choosing a particular dye medium, it is thus necessary to 
check whether it fulfills the Kennard-Stepanov relation. For this it proves helpful to introduce 
a spectral temperature, which is derived using a formal solution of eqn ([T]). Consider the 
Kennard-Stepanov relation for two different frequencies u and u + 6u. From eqn ([T]), we 
obtain 

We can now solve this equation for the frequency dependent spectral temperature Tspec(w): 

h5uj 

Tsv,ec\^) = 7 ^ (23) 



7, a In Bi2{uj) ' 
"'B B2i{u)) 

and eventually perform the limit 5uj — )■ 0. Thus, if the spectral profiles of absorption and 
emission are known, one can easily calculate the corresponding spectral temperature. We 
note that the particular form of Tgpec given by eqn (1231) and (12^ has the advantage of being 
independent of the statistical weights w<^^-^ and the frequency of the zero-phonon line Wq? which 
might not be known in all cases. Moreover, in this form it is sufficient to know the coefficients 
-Bi2,2i(i^) only up to a proportionality factor. 

The Kennard-Stepanov relation holds for a given dye, if its spectral temperature coincides 
with the thermodynamic temperature of the solution, Tgpec (w) = independent of the fre- 
quency uj. Figure H] shows the spectral temperatures, here as function of the wavelength, 
for eight different dye species. They were calculated with eqn. (125]) using the correspond- 
ing spectra of the database of Ref. [42,]. The dyes on the left hand side of the figure fulfill 
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the condition Tspec(w) ~ T to good approximation, while the dyes on the right hand side 
show significant deviations. Noteworthy, the spectral temperature here has a tendency to be 
higher than the ambient temperature, see also 35], |43|. As explanation for the observed de- 



viations, one mainly finds two lines of reasoning in the literature, either based on incomplete 



35 



44j-l46| ■ Our experiment uses 



rovibronic thermalization or on inhomogeneous broadening 
either perylene-dimide or rhodamine 6G dye, for both of which the quantum efficiency is 95% 
or above, and the Kennard-Stepanov relation is well fulfilled. 



III. THERMODYNAMICS OF THE TWO-DIMENSIONAL PHOTON GAS 



Formally, the photon gas confined in the resonator is equivalent to a two-dimensional gas of 
massive bosonic particles. This can be seen by investigating the energy momentum relation. 
We will derive the equilibrium properties of such a gas, and discuss the expected condensate 
fluctuations. 



A. Cavity photon dispersion 

We begin by expressing the energy of a photon in the microcavity as a function of longi- 
tudinal (i.e. along the optical axis) wavenumber kz and transverse wavenumber kr, 

E=-\k\ = -y^k[+kl, (25) 
no no 

where no is the refractive index of the medium. Let us denote the spacing between the two 
cavity mirrors on the optical axis as Dq 1.46 pm in our experiment), and the radii of the 
two curved mirrors, which we assume to be identical, as R (~ Im typically). The boundary 
conditions generated by the mirrors imply a resonance condition 

kzir) = qn/D{r), (26) 



with q as the longitudinal mode number and D(r) = Do — 2(R — ^B? — r^) as the mirror 
spacing at a distance r from the optical axis. In a paraxial limit, with k^ <^ k^ and r <^ R, 
we can approximate the photon energy by the expansion 

£-^ + ^ + i™.,!jV, (27) 

ni 2meff 2 



with nigg = hkz{0)no/c = hcUcutosTT-o/c^ and = (c/no) / y^DoR/2. This yields the dispersion 
relation of a particle with mass nies moving in the transverse resonator that is subject to 



14 



600 
400 
200 


600 
400 
200 


600 
400 
200 


600 
400 
200 




Perylene-diimide 



Phthalocyanine 



520 530 540 
Rhodamine 6G 



510 520 530 
Squarylium 



620 630 640 650 

Coumarin 314 




325 330 335 
Coumarin 343 




435 450 465 
Chlorophyll B 



450 460 
X (nm) 



470 




Figure 4: Spectral temperatures rspec(A) for different dyes at room temperature, derived from tlie 
spectra of the database of Ref. j^]. This numerical determination is most precise in the Stokes region, 
where both absorption and fluorescence are sufficiently strong and thus most precisely known. We 
nevertheless estimate an error of the given spectral temperatures of at least 15 %. The corresponding 
quantum efficiencies are: perylene-diimide: $ = 0.97, rhodamine 6G: ^ = 0.95, squarylium: <1> = 
0.65, coumarin 314: $ = 0.68, phthalocyanine: ^ = 0.60, PPO: <& 1, coumarin 343: ^ = 0.63, 
chlorophyh B: $ = 0.117. 



additional harmonic confinement with trapping frequency Q,. If we in addition account for 
a nonhnear self-interaction of photons that modifies the refractive index following no — )• 
no + n2l{r) , where n2 is a nonlinear index of refraction, and J(r) is the optical intensity, we 
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find that the photon energy is shifted by an amount 

£^int - -mefrc^^/(r) . (28) 

^0 

A non-zero value of n2 can arise e.g. from a Kerr nonlinearity or thermal lensing. Since the 
confined photons have the quadratic dispersion relation of massive particles, we can define 
a thermal de Broghe wavelength by Ath = h/ ^J^wm^kBT, in direct analogy to e.g. a gas of 
atoms. Essentially, Ath is inversely proportional to the average transversal wavenumber sj (k^) 
of the photons, and its value at room temperature is ~ 1.58 pm. Bose- Einstein condensation 
is expected when the phase space density nA^j^ exceeds a value near unity, where n denotes 
the two-dimensional number density. 



B. Statistical theory of two-dimensional Bose gas in trap 

For an exact determination of the phase transition, a statistical multimode treatment is 
necessary. The average photon number at a transversal energy u — h{ui — ^cutoff) is given by 
a Bose- Einstein distribution 

nT,M = , (29) 

e^Bi- - 1 

where g{u) — 2{u/Ml-\-l) is the degeneracy factor (which here increases linearly with energy), 
and the factor 2 accounts for the two possible polarizations. A macroscopic occupation of the 
ground state mode at m = sets in when the particle number N reaches a critical value of 
= Y.u>q'^t,h=q{u). We find 

The typical photon trapping frequency in our setup is ~ 27r ■ 4.1 x 10^° Hz. For room 
temperature, T = 300 K, wc arrive at — 77000, which is experimentally feasible. If 
this critical value is reached, the occupation of transversally excited modes is expected to 
saturate and the ground mode starts to become macroscopically populated. The possibility 
to observe a Bose-Einstein condensation at room temperature can be understood from the 
extremely small effective photon mass rries = fiWcutoHnl/c^, which in our case corresponds to 
2.1eVng/c^ ~ 7 ■ 10~^^kg. This is 10 orders of magnitude below the mass of alkali atoms, 
which extremely enhances the phase transition temperature with respect to atomic systems. 
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Figure 5: Schematic illustration of the excitation exchange between the photon gas and a reservoir 
of electronically excited dye molecules. The photon gas can be seen as an open system, in the sense 
of grandcanonical experimental conditions. 



C. Condensate fluctuations 

In our experiment, photons are initially brought into the system by pumping the dye 
medium with a laser. The pumping is maintained throughout the measurement to compensate 
for photon losses due to unconfined optical modes, finite dye quantum efficiency and mirror 
losses. Despite pumping, the photon gas is expected to both spectrally and spatially relax to 



thermal equi 
being lost 



ibrium, provided that a photon scatteres several time off a dye molecule before 
17|. By optical pumping we generate and uphold a reservoir of electronic excitations 
in the dye medium that can exchange particles with the photon gas. This experimental 
situation is well described by a grandcanonical statistical ensemble, see Fig. [5l in which the 
photon gas acquires the temperature of the dye solution and the photon chemical potential 
is directly related to the excitation level in the medium, as given by eqn (IT^ . 

The particle exchange with the reservoir of excited state dye molecules is of special rel- 
evance for the expected photon statistics and second order coherence of the condensate. It 
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3- 



is well known that there is no ensemble equivalence for Bose-Einstein condensation 
In particular, the particle number fluctuations of the condensate strongly depend on the un- 
derlying statistical ensemble. For a (micro-) canonical ensemble, which is typically realized 
in atomic BEC experiments, the particle number distribution of the ground state changes 
from Bose-Einstein- to Poisson-like when the condensation sets in. This is accompanied by a 
damping of the number fluctuations. On the other hand, the particle number distribution for 
a grandcanonical ensemble remains Bose-Einstein-like and the number fluctuations remain of 
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the order of the average occupation number, Atiq = \/ ('^o) ~ — ^o- This can also be seen 
in the second order coherence of the condensate. The zero-delay autocorrelation function 
(7*^^^(0) = {uQ^nQ — 1)) / (^o)^ is not expected to drop off to g'''^^{0) = 1, as it does for (micro- 
)canonical BECs. Instead, a bunching behavior with g^'^\0) = 2 is expected even below the 

e very interesting to test for such unusually 
large condensate fluctuations. This regime is not observed in Bose-Einstein condensates of 
ultracold atomic gases, but likely to occur for photonic BECs. 



IV. EXPERIMENTS ON PHOTON GAS THERMALIZATION 



A. Experimental setup 

We begin by describing measurements on the thermalization of the photon gas in the dye- 
filled microresonator, which have been carried out at comparatively small photon numbers, 
i.e. far below the critical photon number. A scheme of the experimental setup used in the 
experiments of Refs. 16|, [iT] of our group is shown in Fig. |6l The optical resonator consists 
of two highly reflecting mirrors with reflectivity above 99.997% in the wavelength region 
500 — 590 nm. This reflectivity translates to a finesse of > 100000 for the empty cavity. The 
mirrors are spherically curved with a typical radius of curvature of i? = 1 m. To allow for 
a cavity length in the micrometer regime despite the curvature, one of the mirrors is cut to 
1x1 mm surface size. Both resonator mirrors are glued onto glass substrates. Additional glass 
prisms are attached to the side to allow for a pumping under an angle of 45°, see Fig. |6l The 
typical (effective) cavity length of Dq ^ 1.46 pm is determined from the resonator free spectral 
range, and corresponds to the q = 7 longitudinal mode. The used dyes are either rhodamine 
6G or perylenedimide (PDI), whose absorption and fluorescenc e sp ectra are shown in Fig. [7]. 



The dyes have a high quantum efficiency between 95% and 97% [50|, and spectral temperatures 
close to the thermodynamic temperature of the dye solution. We use organic solvent, e.g. 
methanol and ethylene glycol, and typical dye concentrations of 1.5 x 10~^ Mol/1 for rhodamine 
6G. A careful filtering of the dye solution is necessary before placing it between the mirrors. 
The microcavity is pumped with a laser beam at a wavelength of 532 nm derived from a 
frequency doubled Nd:YAG laser. The light transmitted through one of the cavity mirrors is 
split into two partial beams by a non-polarizing beamsplitter. One beam is imaged onto a 
CCD-camera for spatial analysis, while the other is sent into a spectrometer for a monitoring 
of the emitted spectrum. We use different types of spectrometers, both commercially available 
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Figure 6: Schematic set-up of the microcavity experiment. The dye- filled microresonator is pumped 
at an angle of 45° to the optical axis. The radiation emerging from the microresonator is examined 
both spatially and spectrally. 
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Figure 7: Relative strengths of absorption, BYi{X)lBYf^^ and emission, B^ii}^ j B'^x^ ^ for the dyes 
(a) rhodamine 6G and (b) perylene diimide. 
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Figure 8: (a) Spectral distribution of the radiation emitted from the microcavity far below the 
critical photon number at T = 300 K and T = 365 K (circles). The spectra are in good agreement 
with the Boltzmann-distributed photon energies (lines). For comparison, a T = 300 K Boltzmann- 
distribution is additionally plotted in the bottom graph (dashed line). (Rhodamine 6G in ethylene 
glycol, concentration 5 x 10~'^mol/l, mirror curvatures Ri = R2 = 1 m, cavity order q = 7) (b) 
Normalized spectral distribution of the radiation emitted from the microcavity (connected circles) 
for 4 different cutoff wavelengths Ac — {570, 590, 615, 630} nm, from top to bottom. In addition, 
Boltzmann-distributed photon energies for T = 300 K are plotted (lines). 



and self-built devices. Care has to be taken when the spectrometer makes use of an entrance 
slit. The coupling efficiency to the spectrometer will be different for the various transversally 
excited cavity modes. This problem can be overcome by placing a diffusing plate in front 
of the entrance slit. However, a more advantageous method is to remove the entrance slit 
altogether. This is possible because the microresonator fluorescence can be sufficiently well 
coUimated even without an additional spatial flltering. In Fig. |6l a self-built version of such a 
device is shown. The shown spectrometer uses a rotating grating and a photomultiplier tube. 



B. Spectral and spatial intensity distribution 

In initial experiments, we have tested for the thermalization of the two-dimensional photon 
gas in the dye-fllled resonator jl7|. Figure [8^ shows spectra of light transmitted through 
one of the cavity mirrors for two different temperatures of the resonator setup (top: room 
temperature, T = 300 K, bottom: T = 365 K). The output power for these measurements was 
(50 ± 5 nW), which corresponds to an average photon number in the cavity of A^ph = 60 ± 10. 
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From this, the chemical potential can be determined by numerical solving ^„>o ^r.^tl^) = 
A^ph, with nT^fj,{u) given by eqn 0291) . For the two data sets in Fig. [S^, we obtain fi/k^T = 
-6.76±0.17 (T = 300 K) and /i/zteT = -7.16±0.17 (T = 365 K), respectively. Clearly, both 
measurements are performed far below the phase transition, which for this finite size system 
is expected to occur once the chemical potential becomes comparable to the trap level spacing 
(i.e. /X ^ —HQ). Because of —fi/k-gT ^ 1, the term —1 in the denominator of eqn (12^ can 
be neglected and the distribution becomes Boltzmann-like. As can be seen from the figure, 
besides a derivation near 532 nm from residual pump light, the measured spectra are in good 
agreement with theoretical expectations over a wide spectral range for the two temperatures, 
respectively. We interpret this as evidence for the photon gas to be in thermal equilibrium with 
the dye solution. In another set of measurement, the cutoff wavelength AcutoS = 27rc/wcutofr 
was varied, as can easily be achieved in our setup by a piezo tuning of the cavity length. 
Corresponding spectra are shown in Fig. IHb, along with theory spectra obtained by assuming 
a Boltzmann distribution of the transversally excited modes. Satisfactory agreement between 
theory and experiment is only obtained for the two upper spectra with cutoff wavelengths near 
570 nm and 590 nm. For the two lower spectra with longer cutoff wavelength, no satisfactory 
thermalization is present, which can be attributed to the weak dye reabsorption in this spectral 
regime. These observations illustrate the importance of both emission and reabsorption for 
the thermalization process. 

We have also monitored the spatial distribution of the photon gas at similarly low pumping 
power. As before, the experimental data can be well exp 



ained by assuming a Boltzmann-like 
population of the transversally excited cavity modes |17l |. Furthermore, a spatial concentra- 
tion of light into the center was observed, which can be seen as consequence of minimizing the 
photon energy in the effective trapping potential induced by the mirror curvature. Technical 
applications of this li ght concentration effect could include the collection of diffuse solar light 



to a central spot 



39 



V. EXPERIMENTS ON PHOTON BOSE-EINSTEIN CONDENSATION 

In subsequent experiments, the dye microresonator was operated at higher pump powers, 
allowing for photon numbers near and above the critical regime. To avoid population of 
triplet states and excessive heat deposition, the pump beam in these measurements was 
acousto-optically chopped to 0.5 ps long pulses with 8 ms repetition time. Since the pulses 
are two orders of magnitude longer than the lifetime of the excited dye molecules and four 
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Figure 9: (a) The connected circles give measured spectral intensity distributions for different optical 
pump powers. The legend gives the intracavity optical power, determining the photon number, in 
units of -Pc.exp = (1-55 it 0.6) W. A spectrally sharp BEC peak at the position of the cavity cutoff is 
visible above the critical power. The observed peak width is limited by the spectrometer resolution. 
The inset gives theoretical spectra, (b) Images of the radiation emitted along the cavity axis, below 
(top) and above (bottom) the critical power. In the latter case, a condensate peak is visible in the 
center. 



orders of magnitude longer than the hfetime of the photons (average time between emission 
and reabsorption), the experimental conditions can be considered as quasi-static. 

Typical experimental spectra of the photon gas are shown in Fig. ^ for different optical 
pumping powers At low pumping and correspondingly low intracavity optical power, 

a broad spectral distribution resembling a Boltzmann-like occupation of modes above the 
cavity cutoff is seen. Near the phase transition, the spectral maximum shifts towards the 
cutoff. When exceeding the critical photon number, a spectrally sharp peak at the position 
of the cavity cutoff is visible. The observed spectral width is limited by the resolution of 
the used spectrometer. The experimental results are in good agreement with theoretical 
spectra based on a Bose-Einstein-like population of cavity modes, see the inset of the figure. 
The observed imperfect saturation of the transversally excited modes for even higher photon 
numbers is attributed to an interaction-induced deformation of the effective trapping potential 
j52|. At the phase transition, the optical intracavity power is Pc,exp = (1-55 ± 0.6) W, which 
corresponds to a photon number of (6.3 ±2.4) x 10^. This value is in good agreement with the 
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theoretical prediction of eqn fl3Up . Notably, a similar measurement with PDI gives the same 
result within the experimental uncertainties. Thus, the obtained results seem to be largely 
independent of the used dye. 

Figure [9]3 shows spatial images of the radiation emitted by the dye microcavity (real 
image onto a color CCD camera) both below (top) and above (bottom) the critical power. 
Both pictures show a shift from the yellow spectral regime for the transversally low excited 
cavity modes located near the trap center to the green for transversally higher excited modes 
appearing at the outer trap regions. In the lower image, a bright spot is visible in the center 
with a measured FWHM diameter of (14±2) pm. Within the quoted uncertainties, this agrees 
with the expected diameter of the TEMqo mode, d = 2^^/h\n2Jm^g^i ~ 12.2 pm, yielding clear 
evidence for a macroscopic population of the transversal ground mode. Further, we observe 
that not only the height of the condensate peak increases for higher photon numbers, but also 
its width 16|]. This mode diameter increase suggests a weak repulsive self-interaction mediated 



by the dye solution. The origin of this is most likely thermal lensing, but in principle it could 
also be due to a microscopic Kerr-nonlinearity in the dye medium. Both effects contribute 
to the nonlinear index of refraction of eqn (128|) . if they are treated on a mean- field level. By 
comparing the observed increase of the mode diameter with numerical solutions of the two- 
dimensional Gross-Pitaevski equation, we estimate a dimensionless interaction parameter 53| 
of ^ = (7 ± 3) X 10~^. This is significantly smaller than the values in the range 10~^ . . . 10~^ 



reported for two-dimensional atomic gas experiments 
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and also below the regime in 

which Kosterlitz-Thouless physics can be expected to be relevant Ssj). An indication for the 
latter would be the loss of long-range order. Experimentally, when directing the condensate 
peak through a sheering interferometer we have not seen any signatures for such a spatially 
varying phase. Thus, the observed condensation does not seem to deviate from the BEC 
scenario of an ideal or weakly interacting Bose gas. 

In further experiments, we have investigated the influence of the resonator geometry on the 
observed critical particle number. For that we measured the critical circulating power for vari- 
ous mirror radii and different mirror spacings. From eqn (130|) we expect that the critical power 
Pc = NJloj'^ /2'Kq = {n'^ / 12) {k^TY {u / he) R should increase linearly with the mirror radius R 
and have no dependence on the cavity order q (at least as long as the two-dimensionality of 
photon gas holds). Figure [TUk and the lower graph of Fig. [TUb give the corresponding exper- 
imental results, which confirm the theoretical predictions, both the scaling and the absolute 
values. The upper graph of Fig. [TOb shows the pump power required to achieve criticality as 
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Figure 10: (a) Critical circulating intracavity power as a function of the effective mirror curvature R. 
(Rhodamine 6G in methanol, concentration 1.5 x 10^'^ Mol/1, cavity order q = 7, pulse duration 0.5 jis) 
(b) Pump power at the phase transition (top) and critical intracavity power (bottom), respectively, 
as a function of the cavity order q. (Rhodamine 6G in methanol, concentration 1.5 x 10^^ Mol/1, 
mirror curvatures Ri = R2 = Im, pulse duration 0.5 jis) 



function of the longitudinal mode number q. Here, a decrease is observed for larger mirror 
spacings. This becomes clear by noting that the pump power absorption, which is of order 
1% of the incident light, increases with the thickness of the dye film. Therefore, less pump 
power is necessary to reach a given circulating power. The experimental data can be well 
modeled by assuming that the absorbed pump power at criticality remains constant when 
increasing q (solid line). Note that a quite contrary scaling was observed in early work on 
'thresholdless' microlasers, for which an increase of the pump threshold with mirror spacing 



has been reported 
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As already discussed, for lower intracavity powers the thermalization process is accom- 
panied by a spatial redistribution of photons towards the trap center. This effect is also 
observed for higher pump powers, and we expect it to provide a critical photon density in 

the center of the cavity even for a displaced pumping spot. We note that this effect has been 

I I 

observed in polariton BEC experiments [13], but is unobserved in lasers. For a correspond- 
ing measurement, the pump beam was displaced SOjim away from the trap center. A set of 
spatial intensity profiles for different values of the cavity cutoff wavelength is shown in Fig. 
[TTl By a tuning of the cavity cutoff one can vary the reabsorption probability, which alters 
the degree of thermalization. The lowest shown profile displays results recorded with a cutoff 
wavelength near 610 pm. The low dye reabsorption in this wavelength range inhibits an effi- 
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Figure 11: Intensity distribution of the photon gas along an axis intersecting the trap center for 
different cut-off wavelengths Ac (indicated on the left hand side). The pump beam (dashed line) is 
located outside the trap center and its position as well as its power are kept fixed during the course 
of the measurement. The top curve shows a photon gas at criticality with a ground state population 
of Nq/N < 1% (Rhodamine 6G in methanol, concentration 1.5 x 10~^Mol/l, cavity order q = 7, 
mirror curvatures Ri = R2 = Im., pulse duration 0.5 ps) 

cient thermalization, and spatial accumulation at the trap center. Thus, the maximum of the 
fluorescence remains at the location of the pump spot (shown by a dashed line). If the cutoff 
wavelength is decreased, the reabsorption probability increases and the intensity profiles get 
more symmetrically distributed around the trap center. The top profile shows data recorded 
for a cutoff near 570 nm. Here, the photon gas here seems fully thermalized. The spatial in- 
tensity profile shows a cusp at the trap center, which indicates criticality (condensate fraction 
< 1%). As the pump beam intensity at the trap center is almost neghgible, we conclude that 
the critical photon density is not directly generated by pumping, but is a consequence of the 
spatial relaxation towards the cavity ground state connected to the thermalization. 

VI. CONCLUSIONS 

We have investigated thermodynamic properties of a two-dimensional photon gas in a dye- 
filled optical microcavity. A fluorescence induced thermahzation process leads to a thermal 
photon gas at room temperature with freely adjustable chemical potential. Evidence for a 

Bose-Einstein condensation of photons was obtained from Bose-Einstein distributed photon 
energies including a massively populated ground state mode, the phase transition occurring at 
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the expected critical photon number and exhibiting the predicted dependence on the cavity 
geometry, and a spatial relaxation process leading to a condensation even for a displaced 
pump spot. 

An interesting question is, in what respect does photon Bose-Einstein condensation re- 
late to lasing. In general, the main borderline between a laser and a photonic Bose-Einstein 
condensate is that lasing is a non-equilibrium phenomenon, while Bose-Einstein condensation 
occurs in thermal equilibrium. Thermalization of the photon gas above a cavity low-frequency 
cutoff can be achieved in media fulfilling the Kennard-Stepanov relation, when recapturing 
(spontaneous) emission in one of the many transversal cavity modes, and tuning to a regime 
with strong reabsorption. In such a system, the condition for the emergence of a macro- 
scopically occupied photon mode, which for a laser usually is given by a (small-signal) gain 
condition in an inverted medium, is replaced by a Bose-Einstein criticality condition, as given 
in eqn (I30p for the case of a harmonically trapped two-dimensional photon gas. In essence, 
this is equivalent to the condition of the phase space density exceeding a value near unity. 
The thermalization process will select the cavity ground state to have the biggest statistical 
weight, which even becomes macroscopic when condensation sets in. 

For the future, we plan to investigate condensate properties in more detail. For example, 
an intriguing question is whether the observed self-interaction is sufficient to cause superflu- 
idity. Also, the coherence properties of the photon condensate should be further explored. 
Initial interferometric measurements have already delivered a lower bound for the first order 
coherence length in the centimeter regime |41|. It will be important to also investigate the 
second-order coherence of the photon Bose-Einstein condensate. This is of particular interest 
because of the unusually large intensity fiuctuations that are expected for grandcanonical 
Bose-Einstein condensates. 
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